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Abstract 

We are concerned with the hyperbolic Keller-Segel model with quo- 
rum sensing, a model describing the collective cell movement due to 
chemical signalling with a flux limitation for high cell densities. 

This is a first order quasilinear equation, its flux depends on space 
and time via the solution to an elliptic PDE in which the right hand 
side is the solution to the hyperbolic equation. This model lacks strong 
compactness or contraction properties. Our purpose is to prove the 
existence of an entropy solution obtained, as usual, in passing to the 
limit in a sequence of solutions to the parabolic approximation. 

The method consists in the derivation of a kinetic formulation for 
the weak limit. The specific structure of the limiting kinetic equation 
allows for a 'rigidity theorem' which identifies some property of the 
solution (which might be non-unique) to this kinetic equation. This is 
enough to deduce a posteriori the strong convergence of a subsequence. 
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1 Introduction 



We consider the hyperbolic Keller-Segel model 



' d t u + div (VS(t, y) g(u)) = 0, f>0, j/Gfi 



u (i = 0) = u G L 1 n L 



(O), < u < 1 a.e, 



< 



(1) 



AS + 5 = u in 



k V5 • riQ = on 30. 



Here, the function g(u) is given by 



g(it) = u(l - u) 



therefore we restrict ourselves to solutions satisfying < u(t,x) < 1. The 
problem is posed on Q, it is any bounded domain in Mr, with C 1 boundary, 
and nn(y) is the outward normal to Q at y G 90. One can also take the 
torus = II^ 1 (0,Tj), with Tj > 0, and with periodic boundary conditions 
in O for S; the results and proofs are the same. Notice in particular that 
the normal flux in the equation on u vanishes on d£l and thus the boundary 
is characteristic; therefore we do not need boundary conditions for u (this 
prevents us from investigating questions which rise specific difficulties, see 
[T9] for instance). 

This model represents the density u(t, y) of cells moving with a collective 
chemotactic attraction through the chemical potential S. Their sensitivity 
is limited by the so-called 'quorum sensing' term (1 — u) in g(u). It enters 
a general class of problems in the description of cells movement ( [TTl [23j [T3l 
[Hl[TTJ[T2l[22l[5]). Usually a diffusion term is added to represent the random 
motion of the cells and the above model corresponds to the small viscosity 
limit which has been advocated by several authors, see |TH1 Ell EH H] for 
these aspects. 

This derivation implies that the system JTJ comes with an entropy struc- 
ture as usual (pHEJ). But extra terms enter in this entropy structure because 
of the space dependency of the flux and this leads to a specific difficulties 
(see [HI [H [21 [8] and the references therein). For any C 2 convex function rj 
(the so-called entropy), we have 



d_ 

dt 



77(11) + div(V5 q{u)) + {u-S)[q- grf] (n) < 



(2) 
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where </(£) := g'{&T]'{Q f° r £ e ^- This accounts for the correct jump 
condition on possible discontinuities of u. But due to the dependency of the 
flux on VS(t,y), the above model poses specific difficulties compared to the 
usual theory of quasilinear scalar conservation laws: no a priori compactness 
is known in dimension larger than 1 (no BV bounds or L 1 compactness), 
contraction principle or uniqueness are not known and averaging lemma for 
the kinetic formulation (see below) do not apply (because the transport is 
mostly one dimensional in the direction VS). Even time continuity in L 1 
does not follow from the method we develop in this paper. As a consequence 
we do not know if the full family of solutions to the diffusion approximation 
converges, but only subsequences. All these questions are left open and seem 
difficult. 

Consequently, our proof relies on the weak limit of the diffusion approx- 
imation of that we study through its kinetic formulation. Passing to 
the limit we obtain a weak form of the kinetic formulation of the hyperbolic 
limit. The main ingredient then is to prove a rigidity theorem for the solution 
which implies that the weak limit is a usual entropy solution and that subse- 
quences converge strongly. The kinetic formulation and the main results are 
presented in the next subsection. The diffusion limit is studied in section [31 
and the rigidity theorem is proved in section [U Finally, we analyze the long 
time behavior of solutions in Section [5l Some technical aspects are left in an 
appendix. 

2 Main results 

Our main existence result is the following: 

Theorem 2.1. The system (HJ has a solution u E L°°(IR + x fi), S G 
L°°(R + ;H 1 (p.)) satisfying < u(t,y) < 1, < S(t,y) < 1, and all the 
entropy inequalities J2|) (in the weak sense, with initial data T](uq)). 

Because our method is based on weak limits as mentioned earlier, it is 
more convenient to use the kinetic formulation of JTJ) (see [TH [El EU [7] for 
the theory of kinetic formulations and recent applications). It is a way to 
represent all the inequalities ([2]) in a single equation on the unknown defined 
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on [0,oo) xSlxl, f(t,y,C) = l^ <u ( t , y ), namely 

m(t, y, £) > a bounded measure on [0, T]xS!xl, VT > 0, 

00 

-AS + S = u:= f f(t,y,£)d£ in Q, VS-n n = on dQ. 



(3) 

This is equivalent to J2|), and one can recover (|2j) from (|3j) using that rj(u) = 
J rj' {^)l^ <u ^ t y ^d^ because we can always take r](u) = for u < without 
loss of generality; see also Section [3] for an alternative derivation. 
The outcome of our proof is the following rigidity theorem: 

Theorem 2.2. Consider a weak solution to the kinetic equation 

'!£ + (£- S)g(0% + 9'm v s ■ V y f + R(t, y,0 = % 

m (t,U,0 > a bounded measure on [0, T] xfixR, VT > 0, 

f(0,y,0 = l€<«o( v ), 

00 

-AS + S = u := f f(t,y,£)d£ in Q, VS-n n = on dQ. 



(4) 

which satisfies the properties 

(i) < f(t, y, £) < 1 and / = 1 for £ < 0, / = /or £ > 1, / is nonincreas- 
ing in £, 

(mJ £/iere eziste a constant C > smc/i i/ia£ < C/(l — /) almost every- 
where, 

(Hi) the measure m vanishes for £ < or £ > 1. 

Then, we have f(t,y,£) = If <it(t, y ) and u{t,y) is an entropy solution to JT]). 

The proof of these two results is given in the next sections. The strategy 
in the following : as in [10], we take a parabolic approximation of ([TJ), and 
we intend to pass to the limit in its solution u £ as the viscosity vanishes. 
However, unlike in [10] and as mentioned earlier, the problem lacks a 
priori compactness bounds for u £ when N > 1. Hence, we rather pass to 
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the limit in a kinetic formulation of the approximate problem. The weak 
limit of the sequence f £ = l^ <u ef t y -\ satisfies equation (j3J), with a remainder 
R which can be explicitly computed in terms of / and which satisfies (ii). 
Thus theorem 12.21 implies in turn that u £ converges strongly to u. Let us 
finally mention that an alternative proof for the local existence of strong 
solutions can be carried out, see [4]; however, as stressed by M. Burger, Y. 
Dolak and C. Schmeiser in [1], their strategy does not yield any information 
on the global existence of weak solutions. 

3 The parabolic limit 

In this section, we introduce and study the approximate parabolic system 
with e > 0: 

' d t u £ + div (VS £ (t, y) u £ (l - u £ )) - eAu £ = 0, t > 0, y G fi 

u e (t = 0) = u 6 L 1 nL°°(0), 0<u <la.e, 
< -AS £ + S £ = u £ in n, (5) 
VS £ • n n = on dU, 
Vu £ ■ tiq = on dQ for a.e. t > 0. 

Our goal is to pass to the weak limit in this system but we first state the 
following result 

Proposition 3.1. There exists a unique solution (u £ , S £ ) € Lf oc (0, oo; H 1 (Q))x 
L°°(0, oo; -ff 1 (f2)) of the problem (J5j) and it satisfies the following bounds : for 
all 1 < q < oo, for all T > 0, there exist constants Ci(N,Q,q), C2(N,Q,T) 
such that 

0<u £ (t,y)<l a.e. on [0,oo) x O, (6) 
0<5 e (t,y)<l a.e. on [0,oo)xfl, (7) 

ll'S' e ||L° o (0,oo;VK 2 '9(n)) < Cl; (8) 

^||Vn e || i 2(( 0>T)X Q) + ||5t5 e || i 2(( 0>T ). H i(Q)) < C 2 . (9) 
And for any C 2 convex function r\, we have with the notation in J2|), 

^V(u £ ) + div(V5 e q(u £ )) + (u £ - S £ ) [q - grf] (u £ ) - eAr](u £ ) < 0. (10) 

Proof. Existence and uniqueness of (u £ , S £ ) are easily proved thanks to semi- 
group techniques. The bounds ((6)) follows from the maximum principle be- 
cause and 1 are solutions for all drifts VS £ , the bound j7} also follows 
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from the maximum principle, whereas (JSj) is the regularizing effect for ellip- 
tic equations with smooth coefficients. 

The first bound on Vu £ in ([9]) is obtained by multiplying by u £ the 
evolution equation on u £ . Eventually, differentiating the equation giving S £ 
with respect to t gives 

-Ad t S £ + d t S £ = -div (VS £ u £ (l - u £ )) + eAu £ , 

and the right-hand side is bounded in Lf oc (0, oo; H~ l (£l)) uniformly in e; the 
second bound of (O follows. 



The entropy inequality (flOl) is obtained by multiplication of the evolution 
equation by rj'(u £ ) and using the chain rule. □ 

Next, we pass to the limit in the system J5)). However, because they do 
not provide strong compactness, the bounds on the sequence u £ are insuffi- 
cient to pass to the limit in the nonlinear term 

VS £ u £ (l-u £ ). 

In [10], for N = 1, strong compactness is obtained thanks to uniform BV 
bounds on the sequence u £ ; however, as we have already pointed out, such 
bounds no longer hold when N > 1. Consequently, we pass to the (weak) 
limit in the kinetic formulation for problem ([5]). Our next goal is to present 
this argument. 

We take rj{u) = (u — £)+ in (fTu]l . with and we differentiate (in the 

distributional sense) the inequality obtained with respect to £. This yields 

^ + (£ - S £ )g(0^ + 9'(0VyS* ■ V y f £ - sA y f £ = ^L, (11) 

where m £ (t,y,^) is a nonnegative measure on [0, oo) x H x M. It can be 
written explicitely in terms of u £ , namely 

m £ (t,y,0 := - {d t (u £ - 0+ + div (VS £ l i<u s{g{u £ ) - <?(£))) 
- (u £ - S £ )l^ <u ,g(u £ ) - eA y (u £ - £)+} 
= e\V y u £ (t,y)\ 2 5(Z = u £ (t,y)). 

Notice that < f £ < 1 almost everywhere, and f £ {t,y,£) = when 
£ > 1, / £ (i,y,£) = 1 when £ < 0. Moreover, m £ (t, y, £) = when £ < 
or £ > 1 (in the sense of distributions), and {m £ (t, y, £)} e >o is a family of 
bounded measures on [0, T] x Q x R, VT > 0. 
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Hence, there exists a subsequence, still denoted by e, and functions u = 
u(t, y) G L°°((0, oo) x fi), / = f(t, y, G L°°((0, oo) x Q x R), 5 = S(t, y) G 
L°°(0, oo; W 2,q {Q)), and a nonnegative measure m = m(t,y,^) such that, 
locally in time, 

u £ ^u w*-L°°, 
f^f w*-L°°, 
m £ —^m w — M , 
S £ ^S in L^O.oojW 1 *^)) 

for all p, 1 < p < oo. 

Thus, we can pass to the limit as e — > in equation j5]). All the terms can 
pass to the limit because they are written as weak-strong products except 
9'(£,)^yS £ " ^yf e which yields an extra term. Indeed, we can write 

= div„ {g'm y S £ f) + (u £ - S £ )g'(0f £ - 
In the sense of distributions, as e — ► 0, we have 

divy {g'(0V y S e f) - div y {g'(Z)V y S f) , 
s £ g'(0f £ ^Sg'(0 f. 

But at this stage, we cannot assert that the weak limit of u £ f £ is uf (but it is 
possible to identify it, see (fl~5l) below). Nevertheless, we know that {u e f £ } E >o 
is bounded in L°°; thus, extracting a further subsequence if necessary, there 
exists a function p = p(t, y, £) such that 

u £ f £ p w* - L°°. (12) 

Consequently, 

g'm y S £ ■ V y f £ - sfftVyS ■ V y f + </(£) (p ~ « /) , 

and / is a solution of 

dtf + (C - S)g(£)d 6 f + sKOVyS ■ V y f + g'(0(p - uf) = %m, 

-AS + S = u(t,y) inQ, (13) 
V y S ■ Vn n = on «9Q, 
f(t = 0,y,0 = lf <ti0 ( w ). 
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Moreover, /, u and m inherit the following properties 

< / < 1 a.e., 
/(*,!/, = when£>l, /(i,y,£) = l when £ < 0, 
m(i, x,£) = when £ > 1 or £ < 0, 

[ [ [ m(t, y, £) di dydS, <oo VT > 0. 

And there exists a nonnegative measure z/(t, y, £) such that ^((0, T) xOxR) < 
oo for all T > and 

dsf(t,x,O = -"(t,x,O<0 (14) 

in the sense of distributions. This follows from the fact that 

^/ £ (t,y,£) = -<5(£-n £ (t,y)). 

This means that we have derived the properties (i) and (iii) assumed in 
Theorem 12.21 The remainder term R is here equal to i?(i,y,£) := g'{Q{p ~ 
u /)(Ay>£)- There remains to derive a formula for p which we do now. 

Let if! G Co°([°> oo) xn),ip 2 € Cg°(M) be test functions. 

Then 



p(t,y,Qtpi(t,y) dtdyd£ 

lim/ / u e (t,x)f%t,x,0(p' 2 (0<Pi(t,y)dtdy<% 
£ ^Jo JnxR 

oo 



o Jn 

oo 



e^0 



lim/ / u £ (t,x)(tp2(u £ (t,x)))tpi(t,y) dt dy 



o Jn 

oo 



f°° f d 

lim/ / — (& 2 (!i))f e (t,x,OMt,y)dtd y dZ 
4(^2(£))/(t, x, OMt, y) dt dy d£ 

JQxR «S 

fit, OMt, y) dt d y d£ 

JQxR «? 
Consequently, 
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Next, we integrate this equation on E (t, y are treated as fixed parameters), 
with the boundary conditions f(t,y,£) = and p(t,y,£) = when £ > 1. 
We get 



We show later (see lemma 14"?!]) that this implies the assumption (ii) on 
R in theorem 12.21 

At this stage we have derived the full kinetic formulation for our problem, 
which means that the assumptions of Theorem 12.21 have been obtained in the 
(weak) limit of solutions to the parabolic equation We can turn to its 
proof. 

4 Proof of the rigidity Theorem 12.21 

The technique introduced in [20j is then to compare / and f 2 in order to 
prove that / only takes the values and 1 almost everywhere. Thanks to 
the monotony assumption in (i) (see (fl4l) ). we can then deduce easily that 
there exists u = u(t,y) such that f(t,x,£) = l^ <u u y \. 

Hence, we multiply (j3J) by 2/ and we formally derive an equation for / 2 ; 
the difference f — f 2 satisfies 



jrAf- f) + (£ - S)g(OW ~ f)+g'm y S ■ V y (f - f) + R(l - 2f) = 



We emphasize that this calculation, and the following, seems entirely formal; 
indeed, since / is not smooth, the chain rule 2dtff = dtf 2 for instance, has 
to be justified. Thus, regularizations in (t,y,£) are necessary in order to 
make the argument rigorous. Those are fairly standard (see [SJEHED]), and 
will be detailed in the Appendix. 

It can be seen in the above equation that the key of our method is 
the assumption (ii) on the term R. In the case when R is equal to R = 
g'(p — uf), with p given by (fl~5j) . the inequality in assumption (ii) is proved 
in the following 

Lemma 4.1. For T > 0, set 




(15) 



%n(l-2/). (16) 



C := limsup||u £ || L oo((0,r) x fi). 
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(Notice that C < 1 here). Then, with p given in (fl5l) . we have 

\p(t, y, - u(t, y) f(t, y,0\< Cf(t, y, (1 - f(t, y, £)) 
for a.e. t G (0,T), y G 17, f G R. 

Proof. From (fT5l) . for almost every (i, y, £) G [0, oo) x 17 x R, 



POO 

U(t, y, + y /(t, y, O de' - u(t, y)/(t, y, 

/ d£'f(t,y,0+ / f(t,y, 
Jo J£ 

-f(t,y,0 f(t,y,Z')d£'-f(t,y,t) f(t,y,Ocie 
JO J£ 



(l-/(i,y,OR' 



f(t,y,0 + 



oo 



(l-f(t,y,0). 



Now, remember that f(t,y,£) is decreasing with respect to £ (recall (fl4l) ). 
Therefore /(t, y, £')</(*, 2/, for and l-Z^y^') < l-/(i, y, £) for 

£' < f. And for i G [0,T], y G 0, f(t,y,£) = Oforf > limsup | |u E | | Loo ((0, T) x 

Eventually, we obtain 



< p(t,y,0 " u(t,y)f(t,y,Z) < lim sup | K | | L ~ [/(l - /)] (t,y,£). 



□ 



Now, we integrate ([16]) on f2 x R (notice that for £ < or £ > 1, / = f 2 ) 
We get 



/ (/-/ 2 ){^[K-S)9«)]+A»S 9 '«) 
l(l-2/)|(/-/ 2 ) 



+2/ m(t,y,Odt:f(t,y,Z)dydZ 



< C (f-f 2 )- 



(17) 
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In the above inequality, we have used the facts that / — f 2 > because 
< / < 1, and < S < 1 by the maximum principle. Also, 

|A5| = \S-u\ < 1, 

\g(0\,W(0\ < i vee[o,i], 

/ m(t, y, <%/(i, y, £) dy d£ = - / m(t, y, £)v(t, y, f ) dy d£ < 0. 
JnxR JrixM, 

Consequently, by Gronwall's lemma, we get 

0< / (f-f 2 )(t,y,OdydH<e ct [ (f-f 2 )(t = 0,y,OdydH 

But f(t = 0,y, £) = l ?<no ( y ), and thus (/ - f 2 )(t = 0) = 0. We deduce that 
/(*> Vi = f 2 (t, y, £) for a.e. (t, y, £), and / = or / = 1 almost everywhere. 
Since / is decreasing in £, / = l^ <u n y \, and it is easily checked that in that 
case, p = uf. Hence / is a solution of 

dtf + (£ - S)g(0d i f + V y S ■ Vyfg'(0 = %m, 

and u is an entropy solution of 

d t u + div y (VSg(u)) = 0, t > 0, y G O. 

Remark 4.1. can 6e checked that 

M<u*(t,y) 1 5<«(i,y) ^* - L °° « e ^ n in ^L(( > °°) x 

Hence, we have proved here the following more general result : let {(u n , 5 n )}n>o 
be a sequence such that f n := l^< Un (t,y) satisfies 

dtfn + (£ - S n )g(£)dtf n + g'(0VS n ■ V/„ = %m n + r n , 
— AjSjj + o n = u n , 
< u n < 1, 
/ n (t = 0) = l 5<M o( y ), 

wref/i m n a nonnegative measure and r n — in the sense of distributions. 

Assume that f n — 1 f(t,y,£) w* — L°° and u°(y) — > n°(y) strongly in 
L 1 (0) as n — > oo. 

T/ien f/iere errisis u = u(t,y) £ L°°((0,oo) x fi) swc/i i/iai / = l^< u (t j2/ ) 
and u n —> u in L 1 ((0, T) x n) /or a// T > 0. ,4nd u is an entropy solution 
of CP). 
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5 Long-time behavior 



We wish to mention here a few simple facts on the long-time behavior of 
a solution (u, S) of system {1]). Our motivation comes from the unusual 
complexity of the behavior exhibited in [4j [10] for this limit. From their 
study there appears to be differences between the large time dynamics of 
the parabolic system ([5]) and the hyperbolic system (Q]). The stability of the 
steady states of (pQ) is also discussed in these references. The numerical sim- 
ulations and formal computations presented in these two papers also convey 
a rather good insight of the long time behavior of solutions. Precisely, the 
numerical simulations indicate that for the hyperbolic system (TJ) in dimen- 
sion one, solutions converge to piecewise constant steady states as time goes 
to infinity. In these steady states, regions of vacuum (uqo = below) are 
separated from regions where cells aggregate (uoq = 1 below) by entropic 
shocks. On the contrary, in the parabolic system j5]), when the parameter 
e is large enough (e > |), solutions converge to the constant solution. And 
when the diffusivity parameter e is small, a metastable behavior occurs : 
solutions first get close to the piecewise constant steady states of the hyper- 
bolic system, and then the regions of cell aggregates (called 'plateaus') move 
slowly and at last merge with one another. 

This section aims at proving that any entropy solution to as built 
in Theorem 12.11 converges (in a sense detailed in the next proposition) to a 
steady state solution: 



Our analysis relies on the energy dissipation inherited from the natural 
free energy structure for the chemotaxis systems ([51 1221 [14] ) 



div (VSoo Uoo(l - Uoo)) = in Q 



< 




in 0, 



(18) 



k VSoo • nn = on Oft, 



4/n«(W) 



i\m\\h<p) 



(19) 



2/ n \VS(t,y)\ 2 u(t,y)(l - u(t,y)) dy > 0, 



and consequently. 




n 



\VS(t, y)\ 2 u(t, y)(l — u(t, y)) dy dt —>■ as r — > oo. 



(20) 
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The equality (1191) is proved in [10J when iV = 1, but the dimension does 
not play a significant part here; we reproduce a short proof for the reader 
convenience. First, notice that 

—AS + S = u, -Ad t S + d t S = d t u, 

and thus, multiplying the first equation by dtS and the second by S, after 
integration by parts we obtain 



f ud t S = [ (VS -Vd t S + S d t S) 
Jn Jn 



Notice also that J n uS = \ \S\ \ 2 Hl . Consequently, 

uS = 2 [ Sd t u 
\ Jn 



dt 



-2 f S div (VSg(u)) = 2 f \VS\ 2 g{u). 
Jn Jn 



Therefore, we are led to 

?(*,!/) 



dS(t,y) f du(t,y) ld n QU 2 



And ()I9I) follows. 

Integrating this equality from t = to t = T, we obtain 



r 



o Jn 



NS\ 2 g(u) = \ 



Thus 



(u(T)S(T) - u(t = 0)S(t = 0)) 

[°° [ \VS\ 2 g(u) < hn\ < +oo 
Jo Jn 1 



< -101. 

- 2 1 1 



and (M follows. 



We can now state our main result 

Proposition 5.1. Let (u,S) be a global weak solution of (JTJ) as in Theorem 
HO . Then, for k £ N, fef u fc (i,y) := u(i + fc,y), 5jfc(t,y) := 5(i + fc,y), 
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t > 0, !/ 6 H. T/ien there exists a subsequence (nk)keN such as when k — > oo, 

u nk (t,y)^ Uoo (y) w*-L°°(R + x ft), (21) 
S nk (t,y)^SM inL 2 (0,T;H\Q)), VT > 0, (22) 

/ / lvs^o\u nk (t,y)-u ao (y)\dtdy-^Q, (23) 
Jo Jn 

where Uqo = u 00 (y) G L°°(J7), Sqo = 5*00(2/) G H 2 (Q) are solutions to (fl~8l) 
and < Uoo(y) < 1, < 5oo(y) < 1 and iVSoolitooU - n^) = 0. 

Proof. First step. Weak convergence £>/(«&, Sk). The bounds < Uk, <Sfc < 1, 
and the elliptic regularity 

ll'S'nJ|wi,2((o,T)xf7) + IIVSnJIw^CO.T)*^ < C Vfc G N 

provide us directly with (j2~T1) . (|22l) after extracting subsequences. 

Second step. The limits Uqq and Soo are independent of time. 

First, since the couple (u nk ,S nk ) is a solution of ([TJ), considering a test 
function (p = <p(t,y) in Cg°((0,T) x F), with y>(0,y) = p(T,|/) = for all 
1/6 0, then 

/ / u nk d t (p = - f ( VS nk -Vipg(u nk ). 
Jo Jn Jo Jn 

The right-hand side goes to as k — ► oo according to (|20j) for any test 
function (f. Thus 

/■T r 

Uoo(t, y)d t (f{t, y) dtdy = 



o 

for any test function (p vanishing at t = and t = T. Consequently, Uqo is 
independent oft : Uoo(t,y) = itoo(y)- 

Furthermore, in the weak limit it holds 

- ASoo + = iioo on(0, T) x U, VS'oo • = on dtt. 

and by uniqueness for this problem, we also have = S'oo(y). 

Third step. The limiting equation on u^. 
We now introduce the notations 

A := {(t, y) G (0, T) x ft; VS nk (t, y)^0}, 
B := {(t,y) G (0,T) x ft; u nfc (i,y) -> 0}, 
C := {(t,y) G (0,T) x ft; u„ fc (t,y) 1}. 
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Then A ((0,T) xO\(iUBUC)) = 0, where A is the Lebesgue measure, 
because from (l20l). we deduce that 



T 

2 



VSjk| (t,y) u k (t,y)(l - u k (t,y)) dt dy -> 



/o Jn 

as — ► oo. Hence there exists a subsequence, still denoted by n k , such that 
\VS nk \ 2 (t,y) u nk (t,y)(l - u nk (t,y)) -> a.e. 
It follows from the above strong convergence results that 

[ T [ 1 A |V5„J 2 ^0 = [ T [ 1a|VSoo| 2 , 
io Jn Jo Jn 

/ / lfiWn fc -> = / / lfiUoo, 

io in io in 

/ / l c (l_ Unfe )^0 = / f lo(l-«oo). 

io in io Jn 

Consequently, 

! [ |V^oo|V«oo) = T / |V5oo|^(«oo) = 0, 

io in in 

and |V5oo| 2 g(uoo) = almost everywhere on f2. Thus VS , 00 5('U OC) ) = and 
in particular 

dlVyiVSoo Uoo(l - Uoo)) = 0. 



Fourth step. Proof of (|23ll . 

We have already proved that as k — > oo, 



/ / |V£ n J 2 u nfc (l - u n J -> / / iV^oopttoofl - Uoo) = 0. 
io in io in 

The above convergence results entails that 

T ! \VS n fu 2 nk ^ [ T [ |V5oo| 2 uL. 
o in io in 

And since VS nk VSoo in L 2 ((0, T) x SI), it follows that 



JJ] 
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Writing 

{ u n k ~ u co) = u n k ~ u oo ~ ^ u n k u co + ^ u oo^ 

and using once more the weak convergence of u nk , we obtain 

[ T ( \VS\ 2 ( Unk - Uoo ) 2 -0. 

Jo Jn 

([23]) follows easily, extracting a further subsequence if necessary. 
Fifth step. Kinetic formulation. 

Let / = f(t,y,t;) be the weak limit of l^<-u n (t,y)- Notice that it is not 
obvious that / does not depend on t. Then according to the previous steps 
and to section [31 / satisfies 

dtf + g{m - SooWsf + g'iOVSoc ■ V/ + g'(0(p - uoof) = dt-m, 

where m is a nonnegative measure and p is related to / by equation (|T5j) . 
Moreover, since u nk (t,y) converges to Uoo(y) a-e. on the set {^V-Soo^) ^ 
0}, we deduce that f{t,y,^) = l^ <Uoo ^ a.e. on {VS^ / 0}, and thus 
(p - Uoo/) = and g(£)dtf = on {VS^ ^ 0}. 

□ 

Remark 5.1. In general, stationary states of ([TJ) are not unique, even when 
entropy conditions are required and the mean value on Q is prescribed. Thus, 
it is not obvious that the whole sequence Uk should converge to a stationary 
state Uqo. 



Appendix 

This Appendix is devoted to the rigorous proof of inequality (fl7[) . Regu- 
larizations by convolution are used in order to justify the nonlinear manip- 
ulations which led to equation ([TBI) , as in [2Ql EE]. We focus on the case 
when Q. is an arbitrary bounded domain in Mr, with a C 1 boundary, and 
VS ■ tiq = on dQ; the case when Q = IT^ 1 (0,Tj), and S satisfies periodic 
boundary conditions, is in fact easier, and can be treated in a similar fashion. 

We take 5 1: 6 2 > arbitrary, and tp x E V(R), tp 2 6 V(R N ), tp 3 E V(R), 
with 

< <pi,<Pi,<P3 < 1, 

/ 991 = / (f 2 = / <p 3 = 1, 

Jr Jr n Jr 
Suppc/^ E [-1,0], Suppv?2 E Bi, Slippy E [-1,1]. 
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We set 5 = (81,62), and 



<ps(t,y,0 = 




and for (t, y, f ) £ [0, 00) x R 7 ^ 1 

fs(t,y,£):=f*<ps(t,y,£) = / / / f5(t',y',a')V5(t-t',y-y',(-Z')dt'dy'dti', 

jr Jn jm. 

rag := m* <p$. 

Then /<5 and ms are smooth functions of t, y, £ for all 5 > 0, and < 
/<5 < 1) JT^j > 0. Moreover, is a solution of 

dtfs + (£ - S)g(0dzf s + g'iO^yS ■ V y fs + R*<ps = d^m s + r s , (24) 
and the remainder r,$ is equal to 

r 5 = (H-S)g(0dtf6-[(t;-S)g(0dzf}*<PS 
+g'(0V y S • V y / 5 - [g'(0VyS ■ V„/] * 

We wish to stress that equation (f24l) holds everywhere in (0, 00) X £1$ X R, 
and not in (0, 00) xUxl, where 

n s := { y G 0, d(y,an) > <5}. 

This yields a small difficulty when integrating equation (|21"1) on (0, 00) x 
f&<5 x R, because V<S ■nci s (y) / on d&s even though V5-n^(y) = on 50. 
However, this difficulty can be overcome by using the regularity of S and of 
the boundary <9f2. 

Before writing an equation for fs — f§, let us first prove that rs — > in 
L 1 ((0,T) x O x (-R,R)) for all T,R> 0. In the rest of the appendix, we 
set z = 6 l^ 2 , with z = i, ^ = for 1 < j < N, z N+ i = £, 

Accordingly, we define the differential operators 

d d d 

do = -wr, di = — 1 < i < N, d N+ i = — 
ot dyi d£ 

and we set Q := (0, 00) x $7 x R. 
Then for instance, we have 

= (£ - s)g(0f * dm - [(£ - s)g(0f] * dm 

+ [%((£- S)g(Z))f\*w 

= [ {G(z)-G(z'))f(z')d N+l vs(z-z')dz' 
JQ 

+ [d N+1 G f] * f 5 
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where 

G(z):=(£-S(t,y))g(0, z = (t,y,0- 

Then [d N+1 G f] * <p s converges to d N+1 G f in L p ((0, T) x O x (R, R)) for all 
T,R> 0, And setting ipk(z) = ZkdN+\<p{z) for 1 < k < N, we have 

/ (G(z) - G{z')) f{z')d N+W s(z - z') dz' 
JQ 

= [ [ d k G(rz + (1 - t)z') f(z')^ s (z - z') dz' dr 
JQ JO 

The above integral converges to 

d k G(z)f(z) [ Mz) dz' 

in L 2 ((0,T) x Q x (-R, R)) for all T, R > (recall that 5 is bounded 
in L°°(0,T; W 2 > g (Q)) n iy 1,2 (0, T; thanks to proposition EU). But 

Jrjv+2 ^fc(z') da/ = if k ^ N + 1 and / r jv+2 ipN+i{z') dz' = — 1. Thus 

(£ - S)g(.ZWs - [(£ - 5)5(0^/] * w 

converges to in L 2 ((0, T)xSlx (-JS, R)) for all T, i? > 0. The other term 
can be treated in a similar way, using the bounds on S derived in proposition 

EH 

We now go back to the equation on fs; since fs is smooth in t,y,£, we 
can use the chain rule and write, for (t,y,£) G (0, oo) x !lj x I, 

dt {fs ~ /|) + (£ - S)g(0dt {fs ~ /|) + </(£)V yl S • V y {f s - /f ) + 

+ R * <p s (l - 2f s ) = dzms(l ~ Vs) + r s {l - 2f s ) 

We now integrate the above equation on Q$ x R; notice that since / = for 
£ > 1 and / = 1 for £ < 0, we have fs - /| = for £ < — S or £ > 1 + (5, and 
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similarly, ms, r$ = for £ < — S or £ > 1 + 6. Thus 

d_ 

It 



[ {fs-fl) = [ (fs-fI)[dt((Z-S)g(0)+g'mS} 
-I R*<p s (l-2fs) 

J xR 

+2 / msdsfs + [ r 5 (l - 2f s ) 

- / (fs- fs) ■ n n6 (y) dS(y)d£ 

JdClgXR 

< C I (f s -f$)+C [ |/-/ 2 |*^ 



+ ll r <s(*)llii(^x(-i.2)) +C / |V5-nn,(y)| dS(y) 

Jdn s 

In the above inequality, we have used the fact that d^fs = —v * tpg < 0, 
where v was defined in (TT4l) . together with lemma [4TT1 Moreover, notice that 
since the function x i— ► x — x 2 is concave, by Jensen's inequality, we get 

/ (f-f 2 )*<Ps< [ {fs - f!) • 

And since S belongs to W 2 ' q for all q < oo, VS 1 G C°' Q (f2) for some < a < 1; 
remember that we have assumed that the boundary 30 is at least C . In 
such conditions, it is easily proved that 

\VS-nn s (y)\ dS(y)^0 

dn s 

as 6 — > for almost every i > 0. In the following, we set 



^(t) = lk«(<)lb(nx(-i,2))+C / |V5-n n ,(y)| dS(y), 

and we have proved that Us — ► in Lj oc ([0, oo)) as 5 — > 0. 
Thus we are led to 



dl 



■f {fs-f 2 s)<cf (f s -fl) + u s (t). 
Jn^xR Jn s xR 

Consequently, by Gronwall's lemma, 

(/«(*) - fs{tf) < e ct [ (f s (t = 0) - fs(t = 0) 2 )+ f e°^U s (s)ds. 
n s xR Jn s xR Jo 

(25) 
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There only remains to prove that fg(t = 0) — fg(t = 0) 2 goes to as 
8 — > 0. This is a consequence of the fact that fg(t = 0) strongly converges 
to f(t = 0) = l£<u , but the latter is not obvious since 

poo r 

f s (t = o, y ,o= / fV,i/,?)M-t>y-i/,Z-?)<tfdj/dg. 

Jo JnxR 



We therefore use the same technique as in [21 J : since the proof is strictly 
identical to the one in |21j . we only recall briefly the main arguments. 
Let 

II* ( s 



then d t T h (t) = —fan and thus for (t,y,£) € (0,oo) x Q s x R, 

fs(t = 0, y,£) can be written as 

fs(t = 0,y,0 = l s {y,0 

where 1 1 /<y 1 1 z,°° < C81/82. Passing first to the weak limit as 81,82 ~~ * 
with 81/82 — ► in the above equation entails that the weak limit F of 
/$(t = 0, y, £) satisfies 

F = d ( M + l i<Uo{y) 

for some nonnegative measure M vanishing for large £. This leads to F = 
-k<«o(y) an( ^ = thanks to a lemma in [21] . Then, the above formula for 
fs(t = 0, y, £) is used once again to find the weak limit of fg(t = 0,y,£) 2 . It 
is easily proved that 

f s (t = 0,y,0 2 ~* F 2 = ls<u (y), 

and thus the convergence is strong. 

Consequently, fg — /| converges to in Lj oc ((0, 00) x O x R). Since 
fs-fg^f- f 2 m 2^.(0, 00; L x (n x R)), we deduce that / = or / = 1 
almost everywhere. The rest of the proof, exposed in section [H is therefore 
justified. 
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